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We describe an extension of the hadronic SU(3) non-linear sigma model to include quarks. As a result, we 
obtain an effective model which interpolates between hadronic and quark degrees of freedom. The new parameters 
and the potential for the Polyakov loop (used as the order parameter for deconfinement) are calibrated in order 
to fit lattice QCD data and reproduce the QCD phase diagram. Finally, the equation of state provided by the 
model, combined with gravity through the inclusion of general relativity, is used to make predictions for neutron 
stars. 



1. Introduction 

Dense matter hadronic models, such as the ones 
used to describe neutron stars, work in a deter- 
mined range of energies |1I2I3] . They can, in ad- 
dition, fulfill some symmetries of the underlying 
theory (QCD) such as chiral symmetry |4I5I6| and 
scale invariance |7l8j but they do not include de- 
confinement to quark matter. On the other hand, 
the broadly used bag-model [9] includes quark de- 
grees of freedom but does not include chiral sym- 
metry. Models such as the quark-NJL and quark 
sigma-models [10|ll|12j include that symmetry, 
but do not directly incorporate hadronic degrees 
of freedom or quark confinement. 

Our goal is to construct an effective model that 
contains hadronic and quark degrees of freedom 
present at different densities/temperatures but 
that can also coexist in a mixed phase. This al- 
lows the deconfinement phase transition to be a 
steep first order as well as a smooth crossover and 
cases in between. The last two possibilities, de- 
spite being predicted by lattice QCD, cannot be 
achieved by the usual procedure that connects, by 
hand, two different models with separate equa- 
tions of state for hadronic and quark matter at 
the chemical potential at which the pressure of 



the quark EOS exceeds the hadronic one [T3"] . 

In order to achieve this goal, we employ a single 
model for the hadronic and quark phases. Our 
extension of the hadronic SU(3) non-linear sigma 
model also includes quark degrees of freedom in a 
spirit similar to the PNJL model [14] , in the sense 
that it is a non-linear sigma model that uses the 
Polyakov loop $ as the order parameter for the 
deconfinement. This is a quite natural idea, since 
the Polyakov loop is related to the Z(3) symmetry, 
that is spontaneously broken by the presence of 
quarks. In QCD the Polyakov loop $ is defined 
via $ = |Tr[exp(i J drA^)], where A4 — iAo is 
the temporal component of the SU(3) gauge field. 

2. The Model 

The Lagrangian density of our non-linear sigma 
model becomes: 



L — LKin + Lint + Lself + L$B — U, 
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where besides the kinetic energy term for 
hadrons, quarks, and leptons (included to insure 
charge neutrality) the terms: 

Lint = - J2i TpiilobjiuU + 9i<$>§ + 9i P np) 

+ml]ip u (2) 
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Lseif = -\{mluj 2 + m 2 p p 2 + m 2 cf) 2 ) 
-. 94 +3a,V + ^ + 2 #) 

+ifco(^ 2 + C 2 + <5 2 ) - fci(a 2 + C 2 + 5 2 ) 2 
-fe (4 + T + 3ff2(52 + C 4 ) 
-k 3 (a 2 ~ 6 2 )( - k 4 ln^g^, (3) 

- mlUa + [fi m \h - ^mlu) C, (4) 

represent (in mean field approximation) the in- 
teractions between baryons (and quarks) and vec- 
tor and scalar mesons, the self-interactions of the 
scalar and vector mesons and an explicit chiral 
symmetry breaking term, responsible for produc- 
ing the masses of the pseudo-scalar mesons. The 
Polyakov-loop potential U will be discussed in 
detail below. The underlying flavor symmetry 
of the model is SU(3) and the index % denotes 
the baryon octet and the three light quarks. The 
mesons included are the vector-isoscalars u> and 
(j>, the vector-isovector p, the scalar-isoscalars a 
and f (strange quark- ant iquark state) and the 
scalar-isovector S. The isovector mesons affect 
isospin- asymmetric matter and are consequently 
important for neutron star physics. The mesons 
are treated as classical fields within the mean-field 
approximation [15] . A detailed discussion of the 
purely hadronic part of the Lagrangian, that de- 
scribes very well nuclear saturation properties as 
well as nuclei properties can be found in |16I4I17| . 

Finite-temperature calculations include a heat 
bath of hadronic and quark quasiparticles within 
the grand canonical potential of the system. Fi- 
nite temperature calculations also include a gas of 
free pions and kaons. As they have very low mass, 
they dominate the low density/ high temperature 
regime. All calculations were performed consid- 
ering zero net strangeness except the zero tem- 
perature star matter case. The reason for such 
exception is that the time scale in neutron stars is 
large enough for strangeness not to be conserved. 

The effective masses of the baryons and quarks 
are generated by the scalar mesons except for a 
small explicit mass term Mq (equal to 150 MeV 



for nucleons, 354 MeV for hyperons, 5 MeV for 
up and down quarks and 150 MeV for strange 
quarks) and the term containing the Polyakov 
field $: 

M B = g Ba <J + 9bst 3 S + g B cC 

+M 0B +g B ^ 2 , (5) 

M* = g qa o + g qS T 3 S + g qC ( 

+Mo 5 +.g g *(l-$). (6) 

With the increase of temperature/density, the a 
field (non-strange chiral condensate) decreases its 
value, causing the effective masses of the hadrons 
to decrease towards chiral symmetry restoration. 
The Polyakov loop assumes non-zero values with 
the increase of temperature/density and, due to 
its presence in the baryons effective mass (Eq. [5]) , 
suppresses their propagation. On the other hand, 
the presence of the Polyakov field in the effective 
mass of the quarks, included with a negative sign 
(Eq. [6} , insures that they will not be present at 
low temperatures/densities. 

The behavior of the order parameters of the 
model is shown in Fig. [T]for neutron star matter 
at zero temperature. The calculations for neutron 
stars assume charge neutrality (which is essential 
for their stability) and beta equilibrium. In this 
case, the chiral symmetry restoration, which is a 
crossover for purely hadronic matter, turns into a 
first order phase transition by the influence of the 
strong first-order transition to deconfined matter. 
The fact that the value of the chiral condensate 
increases at a certain chemical potential is due to 
the fact that it is connected to the baryon den- 
sity value that decreases at the phase transition 
(the baryonic number of quarks is one third of 
the hadronic one). The model is consistent in the 
sense that both order parameters are related. 

The coupling constants related to the Polyakov 
loop in the effective mass formulas (Eq. [5] and 
Eq. [6]) are high but still finite. This causes the ef- 
fective masses of the degrees of freedom not effec- 
tively present in each phases to be high but also 
finite (Fig. [2]). The effective normalized masses 
of the particles are directly related to the onset 
these particles appearance in the system. 
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Figure 1. Order parameters for chiral restoration 
and deconfinement to quark matter for star mat- 
ter at zero temperature. 



Figure 2. Effective normalized mass of baryons 
and quarks for star matter at zero temperature. 



The potential U for the Polyakov loop reads: 

U = (a T 4 + a lP L 4 + a 2 T 2 ^ 2 )^ 2 

+a 3 T 4 In (l-6$ 2 + 8$ 3 -3$ 4 ). (7) 

It is based on [18119] and adapted to also include 
terms that depend on the chemical potential. The 
two extra terms (that depend on the chemical po- 
tential) are not unique, but the most simple nat- 
ural choice to reproduce the main features of the 
phase diagram at finite densities. The coupling 
constants for the baryons (already shown in |17| ) 
are chosen to reproduce the vacuum masses of the 
baryons and mesons, nuclear saturation proper- 
ties, and asymmetry energy as well as the hyperon 
potentials. The vacuum expectation values of the 
scalar mesons are constrained by reproducing the 
pion and kaon decay constants. The coupling con- 
stants for the quarks (g qu! — 0, g qt f, = 0, g qp = 0, 
g qiy = -3.00, g qS = 0, g q< = -3.00, T = 270/200 
MeV, a a = -1.85, a x = -1.44xl0~ 3 , a 2 = -0.08, 
a 3 = -0.40, g B $ = 1500 MeV, g q<s> = 500 MeV) 
are chosen to reproduce lattice data as well as 
known information about the phase diagram. 

More specifically, the parameters ao and a 3 of 
the Polyakov potential were fit to reproduce the 
same pressure functional as the one in the model 
presented in Rcfs. |18I19| for the quenched case 



for the range of temperature of interest (Fig. 
In this way we ensure the agreement with lattice 
at the zero chemical potential limit. The param- 
eter To was chosen to be 270 MeV for quenched 
calculations reproducing a very strong first order 
phase order transition at 270 MeV temperature 
(as in Ref. [18]). The parameter To was rescalcd 
for the calculations including quarks and hadrons 
also as in Ref. [IS]- In our case it changed from 
To = 270 to 200 MeV reproducing, at zero chem- 
ical potential, a peak on the susceptibility for the 
Polyakov loop at T = 190 MeV (Fig.©. 

The functional form of the potential in our 
model is very similar to the one in Refs. [T8 19 
with a leading term proportional to temperature 
to the fourth power. At the moment, the other 
powers of temperature (third and second) were 
not considered to avoid a large number of pa- 
rameters. In our case, the potential also con- 
tains a term proportional to chemical potential 
to the fourth power, with parameter a 2 (chosen 
to reproduce a phase transition from hadronic to 
quark matter at a value of four times saturation 
density for star matter at zero temperature), and 
a mixed term between temperature and chemi- 
cal potential, with parameter a x (chosen to re- 
produce a critical endpoint at chemical potential 
fi c = 354 MeV and temperature T c = 167 MeV 
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Figure 3. Pressure for the quenched case. 



for symmetric matter and zero net strangeness in 
accordance with lattice data from Fodor and Katz 
Ref. 20 ]). Together, these parameters lead the 
model to reproduce the QCD phase diagram at 
large densities. More complicated structures for 
the Polyakov potential will be analyzed in future 
work. 

As can be seen in Fig. [5] the transition from 
hadronic to quark matter obtained is a crossover 
for small and vanishing chemical potentials. Be- 
yond the critical end point, a first order transition 
line begins. The critical temperatures for chiral 
symmetry restoration coincide with the ones for 
deconfinement. Since the model is able to repro- 
duce nuclear matter saturation at realistic values 
for the saturation density, nuclear binding energy, 
as well as the compressibility and asymmetry en- 
ergy, we also reproduce results at low densities for 
the nuclear matter liquid-gas phase transition. 

3. Application to Neutron Stars 

In the same way we used lattice QCD data to 
calibrate and test the model for high tempera- 
tures and low chemical potentials, we can cal- 
culate neutron star properties to test the model 
in the high-density /low-temperature regime. It 
is important to notice that up to this point, 
the charge neutrality was considered to be local, 



Figure 4. Order parameters for chiral symmetry 
restoration and deconfinement to quark matter at 
zero chemical potential. 



meaning that each phase was separately charge 
neutral. At finite temperature the two phases 
contain mixtures of hadrons and quarks, which 
are dominated by hadrons or quarks, depending 
on the respective phase. At vanishing tempera- 
ture there is no mixture, i.e. the system exhibits 
a purely hadronic or purely quark phase. 

It is important to notice that the quarks are 
totally suppressed in the hadronic phase but the 
hadrons are suppressed in the quark phase until 
a certain chemical potential (above 1700 MeV for 
T = 0). This threshold, which is higher than the 
density in the center of neutron stars, establishes 
a limit for the applicability of the model. 

For a more realistic approach we allow the 
two phases to be charge neutral when combined 
(global charge neutrality) following 21 . In this 
case, the particle densities change in the coexis- 
tence region causing quarks and baryons to ap- 
pear and vanish in a smoother way (Fig. [6]). 

The density of electrons and muons is signifi- 
cant in the hadronic phase but not in the quark 
phase. The reason for this behavior is that be- 
cause the down and strange quarks are also neg- 
atively charged, electrons are not necessary for 
charge neutrality, and only a small amount of 
leptons remains to ensure beta equilibrium. The 
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Figure 5. Phase diagram. The lines represent first 
order transitions and the circles mark the critical 
end-points. 



Figure 6. Population (baryonic densities) for star 
matter at zero temperature using global charge 
neutrality. 



hyperons, in spite of being included in the calcu- 
lation, are suppressed by the appearance of the 
quark phase. Only a very small amount of A's 
appear right before the phase transition. The 
strange quarks appear after the other quarks and 
also do not make substantial changes in the sys- 
tem. 

The possible neutron star masses and radii are 
calculated by solving the Tolmann-Oppenheimer- 
Volkof equations |22l23j . The solutions for 
hadronic (same model but without quarks) and 
hybrid stars are shown in Fig. where besides 
our equation of state for the core, a separate equa- 
tion of state was used for the crust [51]. The 
maximum mass supported against gravity in our 
model is 2.1M® in the first case and approxi- 
mately 2.0M Q in the second. The predicted radii 
are in the observed range, being practically the 
same for hadronic or hybrid stars. 

Although the stars that contain a phase of pure 
quark matter in our model are not stable, stars 
that contain a core of mixed phase (surrounded 
by a hadronic shell) are. In this case, the mixed 
phase could extend up to approximately 2 km of 
radius but with only a small quark fraction. 



4. Conclusions 

In spite of the fact that our model is relatively 
simple, it is the only one able to take into account 
different degrees of freedom and consequently al- 
low steep as well as smooth transitions between 
different phases. The model is in accordance with 
lattice QCD data and the phase diagram it re- 
produces is able to describe a broad variety of 
regimes, from compact stars to heavy ion colli- 
sions. Calculations along this line are in progress 

Eg. 

We conclude that our model is suitable for 
the description of neutron stars, since it predicts 
masses and radii in agreement with observations. 
Although it does not predict stable stars contain- 
ing pure quark matter, the model allows stars 
that contain a core of mixed phase that can ex- 
tend up to approximately 2 km of radius. Even in 
this case, the reduced maximum star mass is still 
higher than the most massive pulsars observed. 

A major advantage of our work compared to 
other studies of hybrid stars is that we can study 
in detail the way in which chiral symmetry is 
restored and the way deconfincment occurs at 
high temperature/density. Since the properties 
of the physical system, e.g. the density of parti- 
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Figure 7. Mass-radius diagram 



cles in each phase, are directly connected to the 
Polyakov loop it is not surprising that we obtain 
different results in a combined description of the 
degrees of freedom compared to a simple match- 
ing of two separate equations of state. 
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